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ABSTRACT 

In  disc  laser  amplifiers,  the  achievable  inversion  (and  thus  the  gain  and 
energy  output)  is  limited  by  fluorescence  amplification  and  parasitic 
oscillation.  In  this  paper,  analytic  results  and  Monte  Carlo  calculations 
of  both  these  effects  are  presented  for  spheres,  circular  discs,  and 
elliptical  discs.  The  effects  of  gain,  fluorescence  line  profile,  refractive 
indices,  and  disc  thickness  on  fluorescence  amplification  are  presented. 

The  loss  rate  is  seen  to  increase  rapidly  above  an  across-diameter  gain  of 
exp  (5)  or  so,  independent  of  thickness.  Parasitic  oscillation  is  analyzed 
as  a  function  of  edge  and  face  index  ratios,  and  of  edge  roughening.  Above 
exp  (3)  gain,  parasitic  suppression  is  difficult,  even  with  a  rough  edge. 
Parasitics  thus  appear  to  present  more  difficulties  than  fluorescence 
amplification.  The  effects  of  both  processes  on  peak  inversion  achieved 
during  a  pumping  puls^  are  also  presented.  In  conclusion,  the  maximum 
practical  size  of  NdjT  glass  laser  discs  is  estimated  to  be  about  30  cm 
due  to  the  effects  of  parasitic  oscillation  and  fluorescence  amplification. 
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I.  INTRODUCTION 

Multi-disc  laser  amplifiers  are  presently  of  interest  in  the  production 
of  high-energy  light  pulses  for  the  generation  of  dense,  high- temperature 
plasmas  for  x-ray  generation  and  controlled  thermonuclear  reactions.  Disc 
amplifiers,  as  compared  to  rod  amplifiers,  have  the  advantages  of  more 
uniform  pumping,  easy  heat  removal,  small  size  of  individual  glass  pieces, 
and  less  susceptibility  to  self-trapping.  However,  the  design  of  these 
amplifiers  must  ba  carefully  carried  out  if  costly  and  time-consuming  mis¬ 
takes  are  not  to  be  made. 

One  of  the  problems  which  must  be  understood  before  disc  amplifier 

(1-9) 

design  is  undertaken  is  that  of  fluorescence  amplification  (sometimes 

called  amplified  spontaneous  emission,  super fluorescence,  or  -  incorrectly^10^ 
-  superradiance).  Fluorescence  amplification  makes  energy  storage  in 
lasers  more  and  more  difficult  as  the  inversion  level  rises.  Section  II 
of  this  paper  is  a  discussion  of  analytical  and  Monte  Carlo  calculations 
of  the  limitations  due  to  fluorescence  amplification  in  various  geometrical 
shapes  relevant  to  the  disc  laser  problem,  and  of  the  effects  of  varying 
size,  shape,  gain,  refractive  index,  and  line  profile. 

Another  problem  which  arises  in  disc  amplifier  design  is  parasitic 
oscillation.  Such  oscillation  sets  a  sharp  upper  limit  to  the  achievable 
inversion,  and  often  is  a  more  serious  problem  than  fluorescence  ampli¬ 
fication.  In  Section  III  of  this  paper,  parasitic  oscillation  limits  to 
achievable  gain  in  a  single  disc  will  be  calculated  as  a  function  of  the 
reflections  due  to  the  refractive  indices  of  a  disc  and  its  surroundings, 
and  as  a  function  of  the  edge  roughness  of  the  disc. 
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Both  fluorescence  amplification  and  parasitic  oscillation  act  to  limit 
the  maximum  practical  si2e  of  the  discs  in  a  disc  amplifier,  since  both  pro- 
cesse;  get  worse  as  a  disc  becomes  larger  (assuming  constant  inversion  density) . 
Accurate  analysis  of  both  effects  is  therefore  required  to  determine  the  exact 
limitations  on  size  and  gain. 


II.  FLUORESCENCE  AMPLIFICATION 

>?s:nv  materials,  when  properly  pumped,  exhibit  fluorescence  (the  sponta¬ 
neous  emission  of  photons).  The  fluorescence  arises  when  the  upper  level  of 
the  fluorescent  transition  is  populated  by  the  pumping  process,  and 
radiative  transitions  to  the  lower  level  take  place.  Once  the  pumping  is 
turned  off  the  fluorescence  usually  decays  exponentially  with  a  time 
constant  called  the  fluorescent  lifetime  7.  An  exponential  decay  implies 
that  the  process  is  linear  -  that  is,  that  the  fluorescence  is  linearly 
proportional  to  the  population  in  the  upper  level. 

Suppose,  however,  that  the  pumping  was  so  intense  that  an  inversion 
exists;  the  population  in  the  upper  level  is  greater  than  the  population 
in  the  lower  level.  The  pumped  material  then  exhibits  optical  gain  at  the 
transition  frequency,  and  may  be  used  to  amplify  light  in  a  laser  amplifier 
or  oscillator.  However,  the  fluorescence  continues  under  conditions  where 
gain  exists;  in  fact,  the  processes  of  spontaneous  emission  (fluorescence) 
and  stimulated  emission  (gain)  are  inextricably  interconnected.  This  means 
that  fluorescent  light,  once  emitted  within  the  volume  of  the  material, 
will  be  amplified  before  it  reaches  the  edges.  This  increases  the  effective 
fluorescence  loss  rate,  which  is  the  same  as  decreasing  the  instananeous 
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fluorescent  lifetime.  Because  the  gain,  and  thus  the  added  loss,  now  depends 
on  the  inversion  level,  the  fluorescent  decay  process  is  no  longer  linear. 
Instead,  the  decay  rate  becomes  faster  as  the  inversion  level  increases. 

This  fluorescence  amplification  problem  is  especially  severe  in  large  laser 
systems,  since  Che  long  path  iengths  available  lead  to  large  total  gain  and 
thus  large  amounts  of  fluorescence  amplification. 

So  far  we  have  considered  what  happens  when  the  pumping  is  over. 

However,  the  fluorescence  and  gain  processes  also  exist  during  the  pump 
pulse.  If  the  fluorescence  amplification  is  negligible,  so  that  the  system 
is  linear,  then  doubling  the  pumping  amplitude  will  double  the  inversion. 

In  general,  the  peak  inversion  will  be  linearly  proportional  to  the  pump 
pulse  amplitude.  This  proportionality  will  be  destroyed  if  gain  exists  in 
the  material.  In  the  presence  of  fluorescence  amplification,  equal  incre¬ 
ments  of  pumping  will  yield  smaller  and  smaller  inversion  (and  gain) 
increases,  because  the  decay  rate  per  unit  of  inversion  will  increase  with 
the  inversion.  The  details  of  pumping  in  the  presence  of  fluorescence 
amplification  are  studied  in  Section  IV. 

In  this  section  we  will  consider  only  single-pass  fluorescence 
amplification.  That  is,  the  edges  of  the  discs  will  be  assumed  to  be 
totally  absorbing,  so  that  once  the  fluorescent  light  hits  an  edge  it  is 
not  reflected  back  into  the  disc  again.  This  is  the  best  possible  edge 
condition  for  an  amplifier,  since  if  the  disc  edge  is  not  perfectly  black 
the  light  will  be  reflected  back  into  the  disc  and  excite  further  stimulated 
emission,  thus  increasing  the  fluorescence  loss  rate.  We  are  therefore 
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calculating  an  upper  bound  on  gain  due  to  the  effects  of  fluorescence 
amplification,  since  any  actual  disc  will  return  some  reflected  light,  thus 
increasing  the  loss  rate  and  decreasing  the  gain.  We  here  limit  ourselves 
to  the  black-edge  case  because  once  light  is  reflected  from  the  edges 
oscillation  will  occur  at  some  gain  level  and  a  different  method  of  analysis 
must  be  used  to  find  this  level  (see  Section  III). 


A.  Fluorescence  Amplification  in  a  Sphere 

We  first  calculate  the  stimulated  emission  in  a  sphere,  not  because 
sphere  lasers  are  of  interest  but  because  an  exact  analytical  calculation 
is  possible  in  this  case.  In  addition,  the  sphere  forms  a  simple  system  in 
which  to  demonstrate  the  effects  of  line  shape. 

1.  Single-Wavelength  Calculation 

Consider  a  laser  material  in  the  shape  of  a  sphere.  We  will  assume  that 

a  spatially  uniform  inversion  exists  throughout  the  sphere,  and  make  a 

calculation  of  the  instantaneous  rate  of  fluorescence  amplification.  The 

inversion  causes  the  material  to  have  a  gain  coefficient  a,  such  that  the 

intensity  of  a  pencil  of  light  rays  increases  as  P  =  Pq  exp  (a£) ,  where  t 

is  the  path  length  along  the  ray  direction.  Suppose  that  spontaneous  emission 

releases  a  power  of  I  per  unit  volume.  If  the  sphere  has  a  diameter  D,  the 
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total  spontaneous  emission  is  I  n  D  /6.  We  will  calculate  the  amount  by 
which  this  spontaneous  emission  is  amplified. 

Introduce  a  system  of  spherical  coordinates  centered  on  the  surface  of 
the  sphere  (Fig.  1).  The  angle  0  is  measured  from  the  sphere  diameter, 


and  9  gives  the  rotation  around  the  diameter.  The  distance  r  is  measured 
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along  the  direction  defined  by  9  and  0.  The  element  of  volume  in  these 
2 

coordinates  is  r  sin0  dr  d8  d0,  and  the  fluorescence  from  this  volume  is 
2 

I  r  sin0  dr  d9  d0.  We  assume  that  this  fluorescence  is  emitted  isotropically. 
An  element  of  sphere  surface  of  area  dS  at  the  origin  of  coordinates  will 
intercept  a  fraction  of  the  fluorescence  equal  to  the  projected  area  of  dS  in 
the  direction  of  the  volume  element  (dS  OOSQ)  divided  by  the  area  of  a  sphere 
centered  on  the  volume  element  and  passing  through  the  origin  (4  n  r2) . 

However,  the  fluorescence  will  be  amplified  by  a  factor  exp  (etc)  before  it 
reaches  dS.  Thus  dS  will  intercept  a  total  radiation  (spontaneous  plus 
stimulated)  of 


j—  exp  (at)  COS0  sin0  dr  d9  d0  dS. 


Now  integrate  this  expression  over  the  volume  of  the  sphere  to  find  the 
total  flux  on  the  area  dS: 


2r.  it/ 2  DCOS0 


dP 


=  <*9  J r d0  J dr  ear  cas0  sin  0 


x/2  D COS0 


=  I-dS  J~  d0  dr  e0*  COS0  sin0 


o  o 
n/2 


cos 0  sin0  . 
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With  the  change  of  variable  u  =  atf0,  the  flux  becomes 
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2a  e  I  oP  ‘  , 


j  (aD)  J  (aP)' 


By  symmetry,  all  elements  of  the  sphere  surface  receive  the  same  flux,  so  the 
total  emission  hitting  the  sphere  surface  is 
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To  find  the  amount  M  by  which  the  original  spontaneous  emission  has  been 
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multiplied,  divide  by  the  total  spontaneous  emission  1l  x  D  / 6  and  find 


M  =  —  /l-i'W—  -1 

2e  [  f>  V  V  J 


where  we  have  introduced  the  variable  g  =  oP.  Note  that  g  is  the  log  of  the 
gain  straight  across  the  sphere.  The  ratio  of  stimulated  to  spontaneous 
emission  is 


=  a-  [is!  .  i') 
2g  l  p  \  fy 


Our  calculation  has  been  for  a  single  value  of  a.  However,  in  actual 
laser  materials  the  fluorescence  and  gain  are  distributed  in  wavelength 
according  to  some  line  profile,  with  values  from  zero  to  the  peak  value  being 
represented.  The  effective  values  of  the  spontaneous  emission  multiplier  M 
and  the  stimulated-spontaneous  ratio  A  are  therefore  found  by  averaging  the 
single  wavelength  values  over  the  line.  From  the  viewpoint  of  fluorescence 
amplification,  the  worst  possible  case  would  be  a  flat- topped  or  rectangular 
line  profile,  since  in  that  case  the  line  average  is  equal  to  the  peak  value. 
For  more  realistic  line  profiles,  in  which  values  less  than  the  peak  are 
represented,  the  line  average  is  less  than  the  peak  a  by  an  anount  which 
depends  on. the  line  shape. 

It  should  be  noted  that  M  =  1  +  A  for  any  line  shape,  since  by  definition 
M  =  (J  +  K)/J  =  1  +  J/'K  and  A  =  J/K,  where  J  is  the  spontaneous  radiation 
and  K  is  the  stimulated  radiation. 

2.  Lorentzian  Line  Profile 

Let  us,  for  example,  calculate  the  values  of  M  and  A  for  a  Lorentzian 
line  profile.  The  gain  (and  the  fluorescence,  which  has  the  same  line 
shape)  have  the  form 


where  ap  is  the  gain  at  the  peak  of  the  line,  is  the  center  wavelength  and 
W  is  the  full  width  at  half  maximum. 
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Since  all  that  matters  when  averaging  over  a  line  is  the  fraction  of  the 
line  representing  different  values  of  a-  we  need  to  find  a  function  which  gives 
this  amplitude  weighting.  For  any  synsaetric  shape  of  a-  we  nave 


M  = 


1 


r. 

a  M  dX 


/ 


+  00 
a  dX 


M  4^  da 
aa 


da 


Introducing  the  variable  x  = 


we  have 


where  the  Lorentzian  weighting  function  p  (x)  is  given  by 

li 
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PL(x)  = 


dX 

da 


dX 

*  n  - 


r* 


da 


da 


_  Jx(l-x) 


/ 


dx 


yx(l-x) 


.  ■  .-i 


Since  the  integral  in  the  denominator  is  the  beta  function  B(%,%)  =  it, 
we  have 


p  (x)  =  -  1  - 

it  „/x(l-x) 

The  line- averaged  value  of  the  stimulated/ spontaneous  ratio  is  given  by 

\  V  A  (ttpX)  PL(X)  dX 


where  (3  =  OpD-  Unfortunately,  the  various  terms  of  the  integrand  lead  to 
(cancelling)  infinities.  We  therefore  calculate  by  a  less  direct  method. 

Firstly,  we  determine  the  power  series  expansion  of  A  by  expanding  the 
exponential  and  combining  terms.  We  find 


fil  +  !x!  + 

48  280 


where 


(nf3) (n+1) ! 
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Then  the  line  average  is 


(l-x)*  (»3)W)! 


We  reverse  the  order  of  summation  and  integration,  and  use  the  fact  that  the 
integral  over  x  is  a  beta  function  to  find 


.  =  36  .  3e2  ,  5B3  ,  B4  . 
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(n+3)  (ni-1) 


Since  the  series  must  be  summed  numerically  to  find  the  ^  corresponding  to 
some  specific  (3,  the  error  due  to  truncating  at  any  n  must  be  known.  Let 
the  relative  error  due  to  ignoring  terms  beyond  m  be  £  ,  so  that 

^(1,  ^")  K) 

where 

_  3  n  (2r~l 

3n  “  (tri-3)(m-l):  J.‘=1  \  2r 


Then  we  have 


6  = 
m 


(t 

\n=mfl 


As  m  becomes  large,  the  ratio  of  successive  coefficients  of  $  in  the  expansion 
of  A^  becomes 

anfl  _  (nf-3)  C2nH)  1 

ar  (nf4)  (nf  2)  (2nt-2)  +  7  * 

2 

We  therefore  have  an  upper  limit  to  the  truncated  part  of  the  series: 


CD 


E 

n=nrt-l 


<  E 

n=0 


n  4-9/2 


Therefore 


where 


m 

E.”  £, 

n=l 


3.  Gaussian  Line  Profile 

We  may  find  the  effect  of  a  Gaussian  line  profile  having  the  form 
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by  the  same  method  used  for  the  Lorentzian  line.  In  this  case,  the  amplitude 
weighting  function  is 


and  term-by-term  integration  of  the  series  expansion  of  A  yields 


n=l 


In  this  case,  as  in  the  Lorentzian  one,  the  ratio  of  successive  coefficients 
approaches  n  +  7/2  as  n  -*  ®,  and  so  the  truncation  error  limit  given  for  the 
Lorentzian  line  profile  also  applies  to  the  series  expansion  of  A  for  a 
Gaussian  line  profile. 

4.  Summary  of  Analytic  Sphere  Results 

Table  I  gives  numerical  values  of  the  ratio  of  stimulated  to  spontaneous 
radiation  A  as  a  function  of  (3  (the  log  of  the  across-sphere  gain  at  the 
line  peak)  for  flat-topped,  Gaussian,  and  Lorentzian  line  profiles  of  the  gain 
coefficient.  For  small  values  of  8,  we  see  that  A  is  linearly  dependent  on 
8  with  a  coefficient  which  depends  on  the  line  shape.  The  rise  of  A  with  8 
becomes  faster  as  8  approaches  unity,  and  is  very  fast  for  values  of  8  much 
greater  than  one.  We  shall  see  later  that  this  behavior  is  typical  of  any 
shape  of  laser  material.  The  fact  that  decay  rates  rise  sharply  for  8  >  1 
means  that  pumping  becomes  more  and  more  difficult  in  this  region. 
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TABLE  I 


Ratio  of  stimulated  to  spontaneous  radiation  A  in  a  sphere  as  a  function  of 
the  log  of  the  across-sphere  gain  at  the  line  peak  for  flat-topped,  Gaussian, 
and  Lorentzian  line  profiles. 


p 

^Vlat 

agauss 

^LORENT 

.01 

.003760 

.002657 

.001879 

.02 

.007540 

.005326 

.003765 

.05 

.01900 

.01340 

.009470 

.1 

.03852 

.02710 

.01913 

.2 

.07917 

.05543 

.03905 

.5 

.2153 

.1484 

.1040 

1. 

.5000 

.3352 

.2326 

2. 

1.396 

.8790 

.5980 

5. 

12.97 

6.497 

4.153 

10. 

593.6 

206.8 

123.4 

20. 

3.457  x  106 

8.097  x  105 

4.670  x  105 

50. 

6.097  x  1018 

8.778  x  1015 

4.992  x  1017 

100. 

39 

7.984  x  10 

OO 

8.055  x  10 

OO 

4.562  x  10 

5.  Monte  Carlo  Calculation 

The  amount  of  fluorescence  amplification  in  a  sphere  was  also  calculated 
by  means  of  a  Monte  Carlo  optical  power  flow  program.  Since  analytical 
results  are  available  in  the  case  of  a  sphere,  it  may  seem  pointless  to 
duplicate  the  results  with  an  approximate  method.  However,  the  Monte  Carlo 
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program  can  be  used  with  geometric  shapes  of  arbitrary  complexity.  Thus  it 
is  much  more  general  in  its  application  than  analytic  methods,  but  of  course 
any  such  program  should  be  checked  in  a  geometry  where  the  results  are  known. 
This  is  the  reason  for  its  use  in  the  case  of  a  sphere. 

The  Monte  Carlo  program  is  called  ZAP.  It  traces  optical  power  flow  in 
systems  which  are  made  up  of  segments  bounded  by  planes,  spheres,  cylinders, 
cones,  conic  surfaces,  and  so  forth,  or  in  combinations  of  such  segments.  The 
material  in  the  segments  is  specified  by  its  refractive  index  and  absorption 
(or  amplification)  coefficient.  The  system  surfaces  may  have  constant  or 
angle-varying  reflection,  both  specular  and  scattering.  Each  optical  ray 
carries  with  it  up  to  two  hundred  wavelength  intervals.  The  bulk  and 
surface  properties  of  the  system  may  vary  with  wavelength,  and  the  optical 
power  in  each  wavelength  interval  will  change  according  to  the  system 
properties  at  its  particular  wavelength.  Rays  are  started  at  random  from 
specified  sources  (surface  or  volume)  with  given  initial  wavelength 
distributions  of  power,  and  traced  through  the  system  to  find  the  deposition 
or  extraction  of  power  in  all  the  system  elements.  A  ray  cut-off  scheme  is 
used  which  assures  that  ray  truncation  error  is  small  without  wasting  time 
following  weak  rays.  ZAP  was  designed  for  accurate  modelling  of  laser  pump 
cavities  (a  job  it  does  well),  but  it  is  also  well  suited  to  the  calculation 
of  fluorescence  amplification  in  arbitrary  geometries. 

The  sphere  geometry  was  specified  to  ZAP,  and  runs  of  2500  rays  were  made 
for  various  values  of  (3  to  achieve  an  accuracy  of  two  or  three  percent. 

The  rays  were  started  at  random  throughout  the  sphere,  to  correspond  to 
uniform  inversion.  All  rays  hitting  the  sphere  surface  were  totally  absorbed. 
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Runs  were  made  for  a  flat-topped  line  (one  wavelength),  for  a  line  with  a 
Gaussian  profile,  and  for  a  line  with  a  Lorentzian  profile.  The  Gaussian  and 
Lorentzian  lines  were  modelled  by  dividing  a  line  profile  into  equal  intervals 
of  wavelength.  The  gain  coefficient  and  starting  power  in  each  wavelength 
interval  were  then  taken  proportional  to  the  value  of  the  gain  coefficient  at 
the  center  of  the  interval.  It  is  necessary  to  go  a  large  distance  out  on 
the  tail  of  the  Lorentzian  line  in  order  to  get  a  proper  weighting  of  low 
values  of  the  gain  (recall  that  the  Lorentzian  amplitude  weighting  function 
approaches  infinity  as  the  amplitude  approaches  zero) .  The  final  result, 
with  the  above  number  and  spacing  of  wavelengths,  is  about  1%  high  for  a 
sphere  when  the  calculation  is  exact.  To  this  systematic  error,  the 
Monte  Carlo  method  adds  the  previously  mentioned  random  error  of  2  or  3%. 

A  comparison  of  the  analytic  and  Monte  Carlo  results  for  a  sphere  is 
shown  in  Fig.  2.  The  agreement  is  quite  satisfactory,  since  the  exact 
result  lies  within  the  Monte  Carlo  error  limits  in  all  cases. 

B.  Circular  Disc 

We  will  now  consider  the  case  of  a  circular  disc,  which  is  of  considerably 
greater  practical  interest.  We  will  first  perform  an  approximate  analytic 
calculation  which  is  exact  for  infinity  thin  discs  at  low  gain,  and  then 
find  the  results  for  thick  discs  and/or  high  gain  by  the  Monte  Carlo  method. 


Analytic  Approximation:  Thin  Disc,  Low  Gain 


Consider  the  spontaneous  fluorescence  emitted  from  a  small  element  of 


volume  in  a  circular  disc  with  a  totally  absorbing  edge.  Some  of  the  radia¬ 
tion  will  strike  the  edge  directly,  while  the  rest  will  hit  the  faces.  Both 
components  will  be  amplified  as  they  travel.  The  radiation  which  goes  to 
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Fig.  2  -  Analytic  and  Monte  Carlo  values  for  the  ratio  A  of 
stimulated  to  spontaneous  loss  as  a  function  of  the  product 
(3  of  the  gain  coefficient  and  diameter  for  a  sphere.  Curves  are 
shown  for  three  different  line  profiles  of  the  gain  coefficient. 


Fig.  3  -  Side  view  of  paths  of  rays  which  are  trapped  by  total 
internal  reflection  in  a  circular  disc.  The  paths  are  equivalent 
to  straightline  paths  in  an  infinite  cylinder,  formed  by 
multiple  reflections  of  the  disc  in  its  faces. 


Fig.  5  -  Geometry  for  the  calculation  of  the  distance  from 
the  point  of  emission  to  the  disc  edge.  The  formula  is  the 
same  for  0  <  n/2  and  0  >  rt/2,  as  shown. 


We  must  integrate  over  the  disc  volume  (i.e.,  over  r)  and  over  solid 

_2 

angle.  The  fractional  disc  volume  at  radius  r  is  8r  D  dr,  and  the  fraction 
of  solid  angle  at  8,  0  is  (4n)  ^  sin  0  d8  d0.  The  angle  0  varies  from  0 

c 

to  n-6c,  where  9£  is  the  critical  angle,  but  by  symmetry  we  can  take  twice 
the  integral  from  0^  to  n/2.  Likewise  we  may  take  twice  the  integral  from 
0  to  it  for  9.  The  above- critical- angle  integral  for  the  ratio  M  of  total  to 
spontaneous  emission  is  then 
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This  is  somewhat  clarified  by  use  of  the  dimensionless  variable  u  =  2  r/D, 
whence 
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For  the  low  gain  condition  aP  «  1,  we  can  expand  the  exponential  in  its 
pother  series  and  produce  a  simpler  integral: 
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The  first  term  in  the  series  yields  the  simple  result  cos 9£,  but  the  second  term 
is  more  difficult.  Sin  0  cancels  in  numerator  and  denominator,  so  the  0  integral 
is  just  (n/2  -  9c),  leaving 


22 


Mj  =  ccs 9c  +  —  -  0  1 J  duj  d9  u  \^j  1  -  u  sin  9  -  u  oasQ  J  . 

'  0  0 

In  the  9  integral,  the  cos  0  portion  integrates  to  zero,  leaving  twice  the 
complete  elliptic  integral  of  the  second  kind 
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so  that 


M.  =  COS  9  + 
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which  by  the  change  of  variable  q  =  u  becomes 
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thus  yielding  the  final  result 


ml  =  COS 


0 

c  n 


where  G  =  .9159656  •••  is  Catalan* s  constant. 


g.  Radiation  above  the  critical  angle 


Let  us  follow  the  radiation  emitted  at  above  the  critical  angle  for 
total  internal  reflection  by  a  volume  element  at  distance  b  from  the  disc  face 
(Fig.  6).  Starting  at  angle  0  from  the  normal  to  the  face,  the  ray  travels 
a  distance  b  sec  0  before  it  hits  the  face.  The  relative  strength  is  now 
exp  (a  b  sec  0).  A  fraction  T  is  transmitted  through  the  face  and  (for  our 
purposes)  lost,  while  (1-T)  is  reflected  and  crosses  the  disc.  The  path  length 
to  the  opposite  face  is  L  sec  0  (recall  that  L  is  the  disc  thickness) ,  and 
so  the  reflected  ray  is  amplified  by  exp  (a  L  sec  0) .  It  then  again  loses 
a  fraction  T  to  the  outside,  and  (1-T)  recrosses  the  disc,  is  amplified,  and 
is  partially  transmitted  and  partially  reflected.  The  bouncing  process 
continues  until  the  ray  hits  an  absorbing  edge.  Let  the  relative  strength 
at  the  first  encounter  of  the  face  be  P  =  exp  (a  b  sec  0) ,  and  define  the 
gain  on  a  pass  all  the  way  across  the  disc  to  be  K  =  exp  (a  L  sec  0).  Then 
the  ray  strength  is  as  follows  (Fig.  7): 
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P(1-T)K 
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and  so  forth.  The  total  amount  radiated 
3  3 

+  PT(l-T)  K  +  ...  which  sums  to  PT/[1  - 


at  the  start 

at  the  initial  face  intercept 
after  the  first  reflection 
after  crossing  to  the  opposite  face 
after  the  second  reflection 
after  the  second  full  crossing 
after  the  third  reflection 
after  the  third  full  crossing 
is  then  PT  +  PT(1-T)K  +  PT(1-T)2K2 
(l-T)K]  if  we  take  an  infinite  number 
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of  terns.  Of  course,  the  actual  ray  will  terminate  after  a  finite  number  of 
bounces,  but  the  error  in  the  infinite-bounce  approximation  is  small  as  long 
as  the  gain  across  the  disc  times  the  reflectivity  is  small,  in  which  case 
the  ray  strength  is  reduced  to  a  negligible  value  in  just  a  few  bounces.  This 
small-gain  case  is  exactly  the  one  we  are  calculating,  so  we  will  use  the 
infinite-bounce  approximation  above,  which  overestimates  the  actual  result. 

We  must  integrate  over  angle  and  position  to  find  the  overall  average 
gain.  The  fractional  solid  angle  at  angle  0  is  sin  0  d0/2,  and  the  fractional 
volume  is  db/L  (note  that  the  r  and  9  integrals  have  vanished  because  there  is 
no  r  or  9  dependence  in  the  infinite-bounce  approximation).  Then  the  above- 
critical-angle  radiation  contributes  approximately 


0  0 


where  we  have  (by  symmetry)  taken  twice  the  0  integral  from  0  to  9£  rather 
than  the  integral  from  0  to  6£  plus  the  integral  from  x  -  9c  to  n.  For 
ah  «  1  we  expand  the  exponentials  in  the  integrand  and  divide  to  get 
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which  integrates  to 


m2  “  (i  - «  0c)  + » L  (i  -  i) (zhr)  *  - 

While  performing  this  calculation  we  have  made  another  approximation:  we  have 
assumed  that  the  transmission  T  is  independent  of  angle,  and  have  ignored  the 
fact  that  T  depends  on  polarization.  In  actual  fact,  T  decreases  as  0 
approaches  the  critical  angle,  and  more  radiation  returns  to  the  disc  to  be 
amplified.  If  we  use  the  normal-incidence  value  for  T,  this  approximation 
underestimates  the  actual  result,  and  (hopefully)  tends  to  offset  the  over¬ 
estimate  due  to  the  infinite-bounce  approximation.  In  any  case,  both 
approximations  get  better  as  the  disc  gets  thinner,  since  then  even  a  ray 
starting  near  the  edge  takes  many  bounces  to  get  to  that  edge,  and  the 
across-disc  gain  is  low  so  even  if  a  ray  takes  more  bounces  to  exit  than 
we  calculate  (due  to  greater  reflectivity)  the  added  radiation  is  small. 

Y.  Total  Result 

Adding  together  the  results  above  and  below  the  critical  angle,  we  find 
that  for  low  gain  (aD  «  1  and  ab  «  1)  and  thin  discs  we  have 

MEl  +  oD(2-^)(2  +  G)  +  ^  -  l)fe(sFV)  • 

This  result  is  for  one  value  of  a ,  and  thus  corresponds  to  a  flat  line 
profile.  Values  of  m  corresponding  to  other  line  shapes  may  be  calculated 
as  in  the  sphere  case. 


Since  the  low-gain  result  in  any  geometry  will  always  have  the  form 
M  =  1  +  Ca  or  A  =  Ca  (where  C  is  a  constant  depending  on  the  geometry),  the 
flat-topped,  Gaussian,  and  Lorentzian  lines  will  always  have  the  low-gain  ratios 
2:j2-.l  exhibited  for  the  sphere  and  cylinder. 

2.  Monte  Carlo  Calculation 

ZAP  was  used  to  find  the  stimulated- spontaneous  ratio  A  in  circular 
discs.  Again  the  inversion  (and  thus  the  gain  coefficient)  was  assumed 
spatially  uniform.  The  reflection  at  the  face  was  taken  as  the  average  over 
polarization.  In  Fig.  8,  we  see  the  results  for  flat,  Gaussian,  and 
Lorentzian  line  profiles  in  discs  whose  thickness  is  207»  of  their  diameter. 
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Fig.  8  -  Monte  Carlo  values  of  the  ratio  A  of  stimulated  to 
spontaneous  loss  as  a  function  of  the  product  {3  of  the  gain 
coefficient  and  diameter  for  a  circular  disc  of  thickness  equal 
to  20%  of  its  diameter.  The  ratio  of  the  refractive  index 
inside  the  disc  to  that  outside  the  faces  is  1.56.  Curves  are 
shown  for  three  different  line  profiles  of  the  gain  coefficient. 
The  low-gain,  thin-disc  approximation  is  shown  for  comparison. 
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The  refractive  index  inside  the  discs  was  1.56  tines  the  index  outside  and  the 
circular  edge  of  the  disc  was  totally  absorbing.  The  A  ratios  are  shown  as  a 
function  of  g  =  a  D,  where  a  is  the  gain  coefficient  at  the  line  peak  and  D 
is  the  disc  diameter.  Also  shown  is  the  snail-gain,  thin-disc  approximation. 
The  approximation  is  seen  to  be  good  below  pa  .1,  even  for  discs  this  thick. 

The  A  curves  have  the  same  general  shape  as  the  curves  for  a  sphere, 
except  for  a  faster  rise  when  A  is  greater  than  one  (this  faster  rise  is  due 
to  the  longer  path  available  in  the  cylinder  due  to  total  internal  reflection 
from  the  faces). 

The  amount  of  fluorescence  amplification  in  a  disc  depends  on  the  disc 
thickness  and  on  the  ratio  of  refractive  indices  inside  and  outside  the 
sphere,  as  well  as  on  the  gain  coefficient.  As  the  cylinder  index  approaches 
the  index  of  its  surroundings,  less  and  less  of  the  light  inside  it  under¬ 
goes  total  internal  reflection,  and  A  becomes  lower.  For  exact  index  matching, 
no  light  returns  into  the  disc  and  A  has  its  minimum  value.  The  variation  of 
A  with  index  ratio  is  shown  in  Fig.  9  for  various  values  of  g.  The  low- 
gain,  thin-disc  approximation  is  also  shown  in  Fig.  9,  and  is  seen  to  match  the 
Monte  Carlo  results  well  except  at  low  index  ratios,  where  the  angle  for 
total  internal  reflection  increases  and  the  infinite-bounce  approximation 
breaks  down. 

The  Monte  Carlo  results  for  the  variation  of  A  with  relative  disc 
thickness  are  shown  in  Fig.  10.  The  values  for  g  =  .1  are  compared  to  the 
low-gain,  thin-disc  approximation.  The  approximation  is  seen  to  be  good  to 
at  least  80%  disc  thickness.  The  most  interesting  feature  of  these  results 
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LOW-GAIN,  THiN-DISC  APPROXIMATION 
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FRACTIONAL  THICKNESS 

Fig.  10  -  The  fluorescence  amplification  in  a  disc  of  1.56  inside- 
outside  refractive  index  ratio  as  a  function  of  the  ratio  of  the 
disc  thickness  to  its  diameter.  The  increase  is  slow  because  the 
principal  contribution  is  from  rays  undergoing  total  internal 
reflection,  and  this  contribution  is  independent  of  thickness.  The 
low-gain,  thin-disc  approximation  is  compared  with  the  p  =  .1  result. 
Note  broken  vertical  scale. 
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is  the  very  slow  increase  of  fluorescence  amplification  with  disc  thickness. 
This  slow  increase  arises  because  much  of  the  fluorescence  amplification  in  a 
disc  is  due  to  rays  which  are  totally  reflected  at  the  faces  of  the  disc,  and 
the  contribution  of  such  rays  is  independent  of  disc  thickness.  Since  useful 
laser  gain  is  roughly  proportional  to  disc  thickness,  discs  should  be  as 
thick  as  possible  within  the  constraints  of  pumping  uniformity,  thermal 
distortion,  self- trapping,  and  fabrication  difficulty.  There  is  also  an 
obvious  advantage  in  trading  a  thick,  weakly-doped  disc  for  a  thin,  highly- 
doped  disc,  since  fluorescent  lifetime  usually  decreases  with  doping. 

C.  Elliptical  Disc 

Analytical  computation  is  difficult  in  the  case  of  an  elliptical  disc, 
and  unnecessary  in  light  of  the  good  fit  of  the  Monte  Carlo  results  to  the 
sphere  and  circular  disc  cases.  Therefore  only  Monte  Carlo  calculations 
were  made  in  this  case.  The  computations  were  for  discs  whose  major  and  minor 
axes  were  in  the  ratio  1:0.5,  which  is  typical  of  laser  glass  tilted  at 
Brewster's  angle  in  air  if  the  beam  is  round.  The  thickness  was  1/7  of  the 
major  axis,  since  this  corresponds  to  an  actual  disc  in  use  in  a  multidisc 
laser  amplifier  in  our  laboratory.  Once  more,  the  inversion  was  uniform, 
the  inside-outside  refractive  index  ratio  was  1.56,  and  the  curved  edge  was 
totally  absorbing. 

The  ratio  of  stimulated  to  spontaneous  radiation,  A,  is  shown  in  Fig.  11 
as  a  function  of  p  =  a  D,  where  D  is  the  disc  major  axis.  The  results  are 
about  70%  of  those  for  a  circular  disc  of  diameter  equal  to  the  ellipse 
major  axis,  for  low  values  of  g.  As  g  increases,  this  fraction  decreases. 
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Fig.  11  -  Monte  Carlo  values  of  the  ratio  A  of  stimulated  to 
spontaneous  loss  as  a  function  of  the  product  R  of  the  gain 
coefficient  and  major  axis  for  an  elliptical  disc  of  minor  axis 
50%,  and  thickness  14.3%,  of  its  major  axis.  The  inside-outside 
refractive  index  ratio  is  1.56.  Curves  are  shown  for  three 
different  line  profiles  of  the  gain  coefficient. 
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The  radiation  absorbed  by  the  disc  edge  and  that  exiting  from  the  disc 
face  were  found  separately  in  the  elliptical  disc  calculations.  These  values 
are  shown  in  Fig.  12  as  a  function  of  B-  We  see  that  the  greater  part  of  the 
stimulated  radiation  increase  at  high  8  goes  to  the  disc  edge,  rather  than 
to  the  disc  face.  Tn  high  energy  systems,  the  edge  absorber  may  thus  have 
to  withstand  large  loadings. 

D.  Discussion 

We  have  studied  the  problem  of  fluorescence  amplification,  in  which  the 
fluorescence  emitted  by  a  material  with  optical  gain  is  amplified  by  stimulated 
emission  as  it  passes  through  the  material,  thus  increasing  the  natural 
fluorescent  loss  rate.  We  have  defined  the  ratio  of  total  loss  to  fluorescence 
M  and  the  ratio  of  stimulated  to  spontaneous  loss  A. 

Fluorescence  amplification  causes  pumping  to  become  more  and  more  difficult 
as  the  stored  energy  (and  thus  the  gain)  increases  in  a  material.  In  effect, 
the  energy  decays  with  an  instantaneous  fluorescent  lifetime  given  by  t/M, 
where  t  is  the  natural  or  low-gain  lifetime.  This  means  that  as  the  pump 
strength  is  increased,  inversion  will  not  increase  in  proportion.  The  resulting 
pump  efficiency  decrease  is  a  slow  function  of  the  pump  energy,  rather  than 
the  abrupt  upper  limit  due  to  parasitic  oscillation  (Section  III).  The  exact 
details  of  the  process  depend  on  the  geometry  of  the  material  and  the 
fluorescent  line  shape,  but  in  general  pumping  is  quite  difficult  when  the 
gain  across  the  longest  dimension  of  the  material  is  about  exp  (5)  —  150 
or  more  (we  assume  that  the  natural  lifetime  is  not  extremely  long,  so  that 
its  decrease  by  a  factor  of  10  is  a  serious  problem).  The  effect  of  fluore¬ 
scence  amplification  on  pumping  will  be  considered  in  detail  in  Section  IV. 


37 


ELLIPTICAL  DISC 

!  x  0.5  x  0.143 


38 


SL^v.—: 


Fig.  12  -  Amount  of  ratio  M  of  total  to  spontaneous  loss  due 
to  radiation  absorbed  at  edge  (Mg)  and  lost  from  faces  (Mp) 
in  an  elliptical  disc,  as  a  function  of  p.  Most  of  the 
increased  loss  due  to  fluorescence  amplification  is  seen 
to  go  to  the  edge  rather  than  out  the  faces. 
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The  line  shape  of  the  fluorescent  emission  (and  gain)  is  important  because 
a  line  with  broad,  low  wings  has  less  fluorescence  amplification  than  a  more 
rectangular  line,  for  the  same  peak  value  of  the  gain.  For  low  gain,  we 
have  shown  that  the  ratio  A  has  the  relative  values  2:j2:l  for  flat  (rectangular), 
Gaussian,  and  Lorentzian  gain  profiles;  Monte  Carlo  calculations  show  that 
these  relations  are  roughly  maintained  at  higher  gains.  Actual  laser 
materials  often  have  complicated,  multi-peaked  line  profiles,  but  in  principle 
we  could  find  the  single-wavelength  value  A^  for  the  geometry  we  are  using 
and  then  (perhaps  numerically)  integrate  this  with  the  amplitude  weighting 
function  corresponding  to  the  laser  material  we  are  using.  In  practice,  it  is 
undoubtedly  sufficient  to  estimate  where  in  the  flat-Gaussian-Lorentzian 
family  the  actual  line  shape  lies  and  to  interpolate  to  find  A  for  that 
line  shape. 

The  geometry  of  the  amplifying  material  strongly  influences  the  amount  of 
fluorescence  amplification.  In  the  simple  case  of  a  sphere,  A  and  M  depend 
only  on  the  product  of  the  gain  coefficient  a  and  the  disc  diameter  D,  and 
of  course  on  the  line  shape.  When  we  consider  a  circular  disc,  we  find  that 
A  and  M  depend  not  only  on  aP>  but  also  on  the  refractive  index  ratio  between 
the  inside  and  outside  of  the  disc,  and  on  the  ratio  of  thickness  to  diameter 
of  the  disc.  An  index  match  between  inside  and  outside  is  desirable  since  it 
leads  to  the  least  possible  fluorescence  amplification.  However,  in  short- 
pulse  high-power  lasers  the  index-matching  mateiial  may  cause  problems  due  to 
self- focusing  damage,  so  this  remedy  is  not  always  available.  The  fluorescence 
amplification  is  also  least  for  vanishingly  thin  discs,  but  of  course  such 
discs  have  insufficient  gain  along  the  path  of  the  actual  laser  pulse  to  even 
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overcome  their  own  surface  losses.  Since  the  increase  of  fluorescence 
amplification  with  thickness  is  slow,  it  is  not  harmful  to  increase  the  disc 
thickness  until  it  is  limited  by  other  considerations.  For  example,  the 
deposition  of  pump  energy  changes  with  disc  thickness,  but  this  question  is 
beyond  the  scope  of  this  paper.  An  elliptical  disc  has  the  same  type  of  index 
and  thickness  behavior  as  a  circular  disc,  but  the  fluorescence  amplification 
is  somewhat  less  than  in  the  circular  case.  Elliptical  discs  have  the  addition¬ 
al  advantages  of  smaller  size  and  hence  lower  cost  (when  compared  to  circular 
discs  tilted  at  Brewster's  angle  to  a  circular  beam)  and  closer  coupling  to 
the  pump  sources. 

III.  PARASITIC  OSCILLATION 

Parasitic  oscillation  takes  place  when  a  material  with  optical  gain  has 
in  it  (or  through  it)  a  light  path  which  returns  on  itself.  Under  this 
condition,  the  material  will  break  into  oscillation  when  the  gain  is  large 
enough  to  overcome  the  path  losses.  Such  undesired  oscillation  will  make 
it  impossible  to  increase  the  stored  energy  in  the  oscillating  volume,  since 
once  oscillation  begins  all  further  pumping  is  immediately  converted  into 
oscillating  power  and  lost.  If  the  mode  of  oscillation  fills  an  appreciable 
fraction  of  the  laser  volume,  it  thus  sets  a  sharp  upper  limit  to  the  available 
stored  energy,  and  therefore  the  gain.  If  the  mode  gain  does  not  vary  greatly 
with  frequency  across  the  gain  line  profile  (or  if  it  oscillates  so  rapidly 
as  to  effectively  sample  all  values  of  gain  in  the  line),  then  the  place 
where  oscillation  will  begin  is  the  line  peak.  This  means  that  the  maximum 
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gain  coefficient  a  determines  the  oscillation  threshold,  so  that  rectangular 
line  profiles  are  the  best  from  the  view  point  of  parasitic  oscillation, 
since  they  maximize  the  available  laser  output  at  the  onset  of  oscillation. 

Parasitic  oscillation  cannot  be  calculated  by  ordinary  Monte  Carlo 
methods,  since  the  power  in  oscillating  modes  typically  flows  in  precise 
directions  at  any  point  in  space.  Unless  the  Monte  Carlo  algorithm  is 
lucky  enough  to  start  a  ray  at  exactly  the  correct  position  going  in  exactly 
the  correct  direction,  it  will  not  find  the  mode  at  all.  Imagine,  for 
example,  the  simple  case  of  two  reflecting  plane  mirrors  1  cm  in  diameter 
spaced  10  cm  apart,  with  a  laser  material  between  them.  In  the  case  of 
such  large  mirrors,  we  can  ignore  diffraction  and  do  a  simple  ray-tracing 
analysis.  The  mode  occupies  the  whole  volume  between  the  mirrors,  so  the 
Monte  Carlo  method  just  has  to  get  the  angle  right.  In  fact,  however,  there 
is  only  one  exactly  correct  direction  -  the  direction  perpendicular  to  the 
mirrors.  Rays  at  a  small  angle  from  the  correct  direction  will  bounce  between 
the  mirrors  a  number  of  times  before  they  miss  a  mirror  and  exit  from  the 
system.  On  the  average,  a  ray  at  an  angle  £  from  the  system  axis  will 
bounce  N  times,  where  N  is  given  by 


N 


where  d  is  the  mi.ror  diameter  and  K,  is  the  inter-mirror  distance.  If  the 
Monte  Carlo  method  starts  rays  at  random  in  angle,  the  probability  of  getting 
N  or  more  bounces  from  a  ray  is 
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The  rapid  decrease  of  p  as  N  increases  means  that  a  very  large  number  of  rays 
must  be  started  to  get  a  ray  with  many  bounces.  For  our  1  cm  by  10  cm 
resonator,  we  must  start  16,000,000  rays  (on  the  average)  to  get  one  with 
100  or  more  bounces.  Clearly,  random- starting  ray- tracing  methods  are  not 
suitable  for  mode  finding  unless  considerable  modification  is  made  to  them. 
Without  such  modification,  a  Monte  Carlo  algorithm  may  give  what  seem  to  be 
perfectly  reasonable  results  even  in  the  presence  of  an  oscillating  mode 
above  threshold. 

Since  our  Monte  Carlo  program  is  unsuitable  for  the  analysis  of 
parasitic  oscillation,  a  simplified  analytic  method  was  used  instead. 

The  method  uses  a  simple  ray- tracing  viewpoint  which  ignores  diffraction 
effects.  In  addition,  phase  is  ignored  and  only  amplitude  is  considered. 

Both  these  simplifications  are  justified  on  the  basis  that  typical  discs  used 
in  large  lasers  are  very  much  larger  than  a  wavelength,  so  that  modes  have  low 
diffraction  loss  and  close  spacing  in  wavelength. 

We  will  first  consider  lossless  modes  which  rely  on  total  internal 
reflection,  then  modes  which  have  some  loss  and  thus  require  gain  in  the 
laser  material  to  achieve  threshold,  and  finally  the  effects  of  a  rough 
disc  euge  on  the  mode  threshold. 

A.  Lossless  Modes  in  a  Circular  Disc 

It  is  possible  for  a  mode  which  is  totally  lossless  to  exist  in  a 
circular  cylinder  because  all  reflections  are  by  total  internal  reflection. 
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Consider,  for  example,  the  disc  shown  in  Fig.  13.  The  disc  faces  are  in 
contact  with  a  medium  with  refractive  index  n^,  the  disc  itself  has  index 
n^,  and  the  (smooth)  edge  of  the  disc  is  clad  with  an  absorbing  material 
of  index  n^  (we  ignore  any  imaginary  component  of  the  absorber  index). 

If  the  edge  coating  index  n^  is  less  than  n^,  then  total  internal  reflection 
is  possible  off  the  inside  of  the  disc  edge,  and  a  ray  can  bounce  losslessly 
around  the  disc  perimeter  in  a  plane  parallel  to  the  faces  (Fig.  14).  In 
addition,  if  the  outside  index  n^  is  less  than  the  disc  index  n^,  a  ray 
can  bounce  from  face  to  face  of  the  disc  by  total  internal  reflection 
(Fig.  15). 

In  general,  a  ray  may  bounce  off  both  the  faces  and  the  edge,  and  the 

ray  path  will  be  complicated  (Fig.  16).  We  can  describe  such  a  path  in  terms 

of  two  angles:  Y  measures  the  angle  of  the  ray  from  a  disc  diameter  (Fig.  14), 

and  8  measures  the  angle  from  a  normal  to  the  faces  (Fig.  15).  The  angle  of 

incidence  on  the  faces  is  thus  9,  and  the  angle  0  of  incidence  on  the  edge 

obeys  casty  =  casV  sin  9.  Thus  as  8  becomes  smaller,  0  becomes  larger.  For  a 

ray  to  be  lossless,  we  must  have  both  9  >  =  sin  ^(n^/n^)  and 

0  >  0c  =  sin  ^"(n^/n^).  This  dual  condition  sets  a  lower  bound  on  Y,  since 

we  must  have  cas'i  <  COSti  /sin  8  .  Thus  we  must  have 

c  c 


which  implies  that  the  lossless  oscillation  exists  only  outside  a  center 

portion  of  the  disc  of  radius  R  =  D/2  sin  Y  ,  where  Y  is  the  minimum  value  of  Y. 
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Fig.  13  -  The  disc  configuration  considered  in  the  parasitic 
oscillation  calculation.  A  circular  disc  of  index  n2  is 
embedded  in  a  material  of  index  n^.  The  edge  of  the  disc  is 
clad  with  an  absorbing  material  of  index  n3» 


He  conclude  that  if  the  edge  coating  refractive  index  is  less  than  the 
disc  index  lossless  modes  can  exist  in  a  circular  disc.  This  means  that 
as  soon  as  there  is  any  gain  at  all  in  the  material  oscillation  will  commence 
and  further  gain  rise  will  be  inhibited.  The  actual  mode  volume  depends  on 
the  edge/disc  and  outside/disc  refractive  index  ratios.  It  is  clearly 
desirable  to  have  the  coating  index  larger  than  the  disc  index,  since  then 
lossless  modes  do  not  exist. 

B.  Lossy  Modes  in  a  Circular  Disc 

Even  if  there  are  no  lossless  modes  in  a  disc,  it  will  still  oscillate  on 
lossy  paths  if  the  gain  is  high  enough.  We  wish  to  find  the  maximum  gain 
achievable  without  parasitic  oscillation,  so  we  must  find  the  mode  with  the 
least  loss. 

Consider  a  path  which  reflects  off  the  disc  faces  by  total  internal 
reflection,  but  which  has  only  partial  reflection  at  the  edge.  We  will  use 
the  same  angles  Y,  9  as  we  used  in  the  previous  section  to  describe  the  ray. 
Let  us  assume  that  the  edge  reflection  is  given  by  R^,  the  larger  of  Fresnel's 
reflection  coefficients.  Then 


where  r  =  sin  ^  (n^/n^  sin  0)  and  0,  as  before,  is  given  by  0  =  cos  ^(n^/n^ 

cos Y)  when  9  is  equal  to  8^.  The  path  length  between  edge  bounces  is 

D  cosV sin  9  =  D  n„ /n,  cos Y,  where  D  is  the  disc  diameter.  To  oscillate,  the 
c  l  1  7 

net  gain  must  be  equal  to  unity,  or 


a  Dn 

-  9  - -  aty  Y 

j~  sin(0-r)  nl 

|^sin(0fr)J  e 


where  a  is  the  gain  coefficient  at  the  line  peak,  0  =  ay-1  (nj/n  COS' *')  and 
r  ~  sin  ^n2^n3  0^*  We  could  now  vary  V  in  this  expression  to  find  the 

minimum  &  required  for  oscillation,  but  it  is  sufficient  to  examine  the 
behavior  at  Y  =  0  and  Y  =  rt/2.  Solving  for  aP,  we  have 


8(y)  =  a  D  = 


2n^fir, 


sin(0t-r) 


sin(0-r) 


1 


n2  COS  Y 


and 


8(0)  =  2  p  im 


Vp2+P2-1  + 

Vc2+p2-l  -  p 


inhere  p  =  and  <j  =  n3/n2-’  while 


,  2 

g(*/2)  =  -■=£==  . 

Vc  -  1 

Vie  have  again  used  the  variable  8  =  aJ). 

The  calculation  is  only  valid  for  ^  <  n2  and  Oj  >  n2,  or  p  <  1  and 
o  >  1.  In  this  region  8(0)  is  always  less  than  $(ji/2),  So  we  find  that  the 
lowest  threshold  modes  are  those  in  which  the  ray  proceeds  across  a  diameter 
of  the  disc,  bouncing  from  face  to  face  at  the  maximum  angle  which  will 
produce  total  internal  reflection.  Such  modes  obviously  fill  the  entire 
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disc  volume.  The  threshold  condition  may  be  written 


Since  the  worst-case  rays  are  those  that  travel  diametrically  across  a 
circular  disc,  we  expect  the  threshold  for  an  elliptical  disc  of  major  axis 
equal  to  the  disc  diameter  to  be  essentially  the  same  as  the  circular  threshold, 
since  the  same  path  is  available.  The  surface  curvature  at  the  reflection 
points  is  of  course  larger  in  the  elliptical  case-  but  this  should  raise  the 
threshold  only  slightly,  since  all  it  does  is  make  the  resonator  somewhat 
more  unstable. 

In  addition  to  the  modes  involving  reflection  off  the  edges,  a  relatively 
thick  disc  may  also  oscillate  between  its  faces.  The  threshold  for  this  process 
is  easily  calculated,  since  the  face  reflection  at  normal  incidence  is 


and  the  gain  from  face  to  face  is  just  exp  (aL),  where  L  is  the  disc  thickness. 
We  must  have  R  exp(ctL)  =  1,  or 


Of  course,  this  result  also  applies  to  elliptical  discs. 
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C-  Combined  Oscillation  Diagram 

We  may  show  the  effects  of  both  lossless  and  lossy  modes  on  a  single 
diagram  in  the  p,  a  plane  (recall  p  =  n^/^  an<*  0  =  n3^n2^*  ^en  lossless 
modes  fill  the  disc  from  a  radius  R  outward,  we  use 


rr  =  i  V1  '  °2 


to  derive 


[■  -  (?)1 


2  2  . 
p  +  a  =  1 


Thus  the  locus  of  constant  R  in  the  p,  o  plane  is  an  ellipse.  This  result 
applies  if  p  <  1  and  a  <  1;  if  p  >  1  then  there  is  no  total  internal  reflection 
off  the  faces  and  the  answer  becomes  independent  of  p,  since  the  oscillations 
lie  in  planes  parallel  to  the  disc  faces  and  ignore  the  condition  at  the 


faces. 


Lossy  modes  obey 


3  =  2p  071 


I  2 ,  2  . 

fc  +p  -1  +  p 

rr~r  7 

'o  +p  -1  -  p 


so  we  may  draw  loci  of  constant  {3  in  the  p,  a  plane  to  show  how  much  gain 
may  be  achieved  before  oscillation  commences.  The  lossy  mode  calculation  is 


,  »v  .i>t  <  /-’ll  ii 


valid  if  g  >  1  and  p  >  1,  and  again  if  p  >  1  the  expression  becomes  independent 
of  p. 

The  diagram  showing  lines  of  constant  oscillation- free  radius  R  and  lines 
of  constant  across-disc  gain  log  f3  is  shown  in  Fig.  17.  From  the  view  point 
of  oscillation  suppression,  the  best  place  to  operate  is  with  both  edge  and 
face  index-matched,  since  oscillation  is  impossible  under  these  circumstances. 
If  other  considerations  make  such  operation  undesirable,  the  degree  of 
oscillation  difficulty  encountered  with  other  choices  of  refractive  indices 
is  easily  evaluated  from  the  diagram. 

D.  Parasitic  Oscillation  with  a  Rough  Edge 

Since  the  oscillatory  modes  considered  above  rely  on  rays  which  follow 
specific  paths,  one  method  of  reducing  oscillation  tendency  would  seem  to 
be  the  roughening  of  the  disc  edge,  so  that  an  initially  parallel  beam  would 
be  smeared  out  in  angle.  This  ought  to  remove  much  of  the  energy  from  the 
mode,  thus  increasing  the  loss.  However,  energy  will  also  be  scattered  into 
the  mode  from  rays  initially  not  involved.  Clearly,  a  quantitative  analysis 
of  this  situation  is  required.  Again  a  Monte  Carlo  approach  is  inapplicable, 
since  a  ray  which  is  randomly  scattered  at  the  edge  has  the  same  low 
probability  of  getting  into  the  high-gain  region  of  the  mode  as  was  the 
case  with  specular  modes.  Of  course,  we  expect  the  active  mode  volume  in 
space  and  angle  to  be  larger  with  a  rough  edge  than  in  the  specular  case, 
but  the  amount  by  which  this  helps  the  Monte  Carlo  method  is  hard  to  estimate. 

The  rough-edge  oscillation  was  therefore  analyzed  by  a  method  related 
to  the  well-known  method  of  analyzing  modes  in  an  optical  resonator.  An 

52 


-  j±  -  -  j&s&MSiz sgattadas. 


i0 


OSCILLATION 
THRESHOLD 
GAIN  LOG 

\ 

- yS=3 


LOSSY 


. _  LOSSLESS 

0.9 

RADIUS 

s - [FRACTION 

0-75|FREE  OF 
LOSSLESS 
OSCILLATION 


LOSSLESS 
OSCILLATIONS 
FILL  ENTIRE 
DISC 


1.0  1.2 


—  --P 

nr.  r 


Fig.  17  -  Combined  diagram  of  lossless  and  lossy  parasitic 
oscillation  as  a  function  of  the  index  ratios  of  a  circular 
disc.  In  the  lower  part  of  the  diagram  (nj/ny  <  1),  the 
fraction  of  the  disc  free  of  lossless  oscillations  is  shown. 
In  the  upper  part  (n^/^  >1),  the  p  required  to  commence 
oscillation  along  lossy  paths  is  shown.  See  text. 


initial  angular  distribution  of  optical  power  was  assumed  at  the  disc  edge, 
and  the  power  was  then  transported  across  the  disc  (with  gain),  reflected 
with  loss,  and  scattered  in  angle.  In  general,  the  distribution  after  this 
process  was  different,  in  both  amplitude  and  shape,  from  the  original 
distribution.  The  gain  coefficient  was  adjusted  after  each  iteration  to 
keep  the  amplitude  constant,  and  the  process  was  reiterated  until  the 
angular  distribution  of  power  had  reached  an  equilibrium  state.  After 
mary  such  iterations,  the  gain  had  the  value  appropriate  to  the  oscillation 
threshold,  and  the  angular  power  distribution  was  that  of  the  oscillating 
mode.  Once  again,  diffraction  and  phase  were  ignored. 

Since  radiation  which  is  totally  internally  reflected  from  the  disc 
faces  has  much  higher  gain  travelling  across  the  disc  than  radiation  which 
is  only  partially  reflected,  the  power  was  assumed  totally  lost  when  it 
was  partially  reflected.  With  this  assumption,  all  points  on  the  disc 
edge  are  equivalent,  and  only  the  angle  of  the  optical  power  is  important 
(as  opposed  to  its  position).  The  hemisphere  of  possible  angles  was 
divided  into  sectors  of  equal  solid  angle,  and  the  resulting  matrix  of 
powers  was  iterated  as  above.  The  edge  roughening  was  simulated  by  taking 
the  unit  vector  along  the  direction  of  the  reflected  ray,  and  adding  to  it  a 
random  vector  with  uniform  distribution  throughout  a  sphere  of  radius 
F  casB,  where  F  2  1  and  0  is  the  angle  from  the  reflected  ray  to  the  normal. 
The  resulting  vector  was  then  renormalized  to  unity.  Various  amounts  of 
surface  roughening  were  simulated  by  varying  F  from  zero  (specular  reflection) 
to  one  (very  rough).  The  edge  reflectivity  R  was  independent  of  angle. 
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The  increase  in  threshold  due  to  increasing  F  is  shown  in  Fig.  18 
for  various  values  of  the  edge  reflectivity.  There  is  obviously  improvement 
over  the  specular  case,  but  it  is  not  as  great  as  might  be  expected.  From 
a  practical  viewpoint,  it  is  very  hard  to  reduce  the  reflectivity  of  a  rough 
surface  to  the  low  values  achievable  by  index-matching  a  smooth  surface, 
because  of  the  many  pits  and  cracks  in  a  rough  edge.  Therefore  one  should 
be  cautious  when  attempting  laser  improvements  by  edge  roughening. 

E.  Discussion 

We  have  shown  that  lossless  modes  exist  for  certain  refractive  index 
ratios  between  the  inside,  outside,  and  edge  of  a  circular  disc.  Even  if 
these  modes  are  suppressed  by  proper  choice  of  the  indices,  lossy  modes 
will  still  exist.  Such  modes  place  an  upper  limit  on  gain,  since  once  the 
gain  reaches  the  mode  threshold  for  oscillation,  no  further  gain  increase  is 
possible  in  the  mode  volume.  Highest  gain  is  achieved  for  close  index 
matching  between  the  disc  and  the  material  in  contact  with  its  faces,  and 
between  the  disc  and  the  absorbing  coating  on  its  edge.  We  have  shown 
what  gain  is  possible  if  proper  matching  is  not  feasible  for  reasons  not 
connected  with  oscillation. 

Even  if  the  absorbing  material  on  the  edge  is  i  good  index  match  to  the 
disc,  there  may  be  areas  of  less- than-per feet  adhesion,  or  chips  and  cracks 
in  the  disc  edge.  Such  flaws  will  have  reflection  larger  than  the  coating- 
disc  interface,  and  in  practice  will  make  the  achievement  of  very  high 
absorption  at  the  edge  difficult.  If  the  higher-reflectivity  flaws  are  small, 
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Fig.  18  -  The  moderate  improvement  in  oscillation  threshold 
achieved  by  various  amounts  of  edge  roughening.  F  is  proportional 
to  the  amount  by  which  the  reflected  rays  are  smeared  out  in 
angle  about  the  specular  direction. 
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diffraction  will  spread  the  reflected  energy,  hut  we  have  seen  that  such 
spreading  causes  only  small  threshold  improvements.  Thus  a  well-index- 
matched  material  must  be  applied  extremely  well  to  a  flawless  surface  if 
a  high  parasitic  oscillation  threshold  is  to  be  achieved. 


IV.  EFFECTS  ON  PUMPING 

Our  discussion  so  far  has  been  confined  to  an  instantaneous  snapshot 
of  the  fluorescence  amplification  and  parasitic  oscillation  processes. 
Lasers  are  usually  pumped  by  relatively  slow  pump  pulses,  and  we  are 
interested  in  knowing  what  decrease  in  efficiency  is  caused  by  the  loss 
processes  we  have  calculated  during  such  a  pulse.  We  will  first  find  the 
effect  of  fluorescence  amplification,  and  then  the  effect  of  parasitic 
oscillation. 

A.  Fluorescence  Amplification 

In  the  low-gain  case,  when  fluorescence  amplification  is  negligible, 
the  stored  energy  decays  with  the  fluorescent  lifetime  t.  The  gain  has  the 
same  behavior,  and  so  we  may  write  a  differential  equation  for  the  across- 
disc  gain  log  g  in  the  form 


B  =  P(t)  -  B/t 


where  p(t)  is  the  pump  pulse  (in  appropriate  units).  At  higher  gain, 
fluorescence  amplification  increases  the  loss  by  a  factor  M,  and  so  we  have 


B  =  p(t)  -  p  M  (B)/t. 


If  we  normalize  to  t  =  1,  and  assume  a  half-sine  approximation  to  the  pump 
pulse,  we  get 


e  = 


up 

Ko 

2T 


-  ei(3) 


0  <  t  <  T 


where  gQ  is  the  area  of  the  pump  pulse  and  T  i.s  the  base  width  of  the  pump 
pulse  (in  units  of  the  fluorescent  lifetime).  In  the  range  0  £  5,  we 

may  make  a  good  approximation  to  M  for  any  of  the  geometries  and  lineshapes 
studied  by  using  the  form 


M  exp  B  +  G  62  +  H  g3 

Note  that  F  may  be  found  from  the  low-gain  value  of  A,  since  A  —  F  g  when 
6  «  1.  We  may  then  find  G  and  H  by  fitting  a  straight  line  to 
(072  M/p  -  F)/g  =.  G  +  Hg. 

Once  M  is  approximated,  we  may  integrate  the  differential  equation  for 
g  until  it  reaches  its  peak  value.  A  plot  of  such  maximum  values  is  shown 
in  Fig.  19  as  a  function  of  the  pumping  gQ,  for  various  values  of  the  pump 
pulse  width.  We  see  that  fluorescence  amplification  is  more  serious  for 
slow,  large-area  pump  pulses  while  faster  or  smaller  pulses  suffer  less 
loss.  In  the  design  of  a  flashlamp-pumped  laser,  the  effects  of  current 
density  on  the  pumping  efficiency  of  flashlamps,  and  the  effects  of  energy 
input  on  flashlamp  life,  should  be  integrated  with  this  decay  information 
in  order  to  choose  the  optimum  operating  point. 

Spatially  uniform  gain  has  been  assumed  in  all  our  calculations  of 
fluorescence  amplification.  However,  the  removal  of  energy  due  to 
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Fig.  19  -  The  maximum  gain  log  (3  produced  by  half- sine  pump 
pulses  of  various  durations  as  a  function  of  the  pulse  area 
Since  the  inversion  is  proportional  to  g,  this  shows  the 
achievable  inversion  for  various  ratios  of  pump  width  to 
fluorescent  lifetime  as  a  function  of  effective  pump  energy 
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fluorescence  amplification  will  be  spatially  non-uniform  in  general.  Thus  the 
gain  will  become  non-uniform  during  a  pumping  pulse  even  if  the  pumping  is 
uniform.  The  analysis  of  this  problem  is  beyond  the  scope  of  this  paper,  but 
it  is  clear  that  proper  tailoring  of  the  spatial  distribution  of  the  pump 
energy  can  reduce  or  eliminate  the  difficulty. 

B.  Parasitic  Oscillation 

Parasitic  oscillation  causes  a  rapid  limitation  of  stored  energy,  rather 
than  the  slow  limitation  we  have  seen  above  for  fluorescence  amplification. 
This  is  because  a  mode  below  its  oscillation  threshold  contributes  little 
to  the  loss  rate,  but  once  the  mode  is  over  threshold  the  loss  rate  is  equal 
to  the  energy  input  rate.  Thus  the  gain  is  sharply  limited  at  the  oscillation 
threshold.  Figure  20  shows  again  the  curve  of  gain  versus  pump  area  for  a 
pump  pulse  of  base  width  equal  to  the  fluorescent  lifetime,  and  also  shows 
the  type  of  gain  saturation  expected  from  the  onset  of  oscillation.  The 
oscillation  threshold  is  shown  for  various  ratios  of  the  refractive  index 
of  the  disc  edge  coating  to  the  refractive  index  of  the  disc;  this  ratio 
determines  the  reflectivity  at  the  edge.  We  have  again  assumed  a  disc-to- 
environment  index  ratio  of  1.56.  This  diagram  illustrates  that  unless 
very  good  index  matching  (or  other  reflection  suppression)  is  used,  parasitic 
oscillation  will  limit  available  energy  storage  well  before  fluorescence 
amplification  losses  become  appreciable. 


I  V.  CONCLUSIONS 


The  effects  of  fluorescence  amplification  and  parasitic  oscillation  on 
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Fig.  20  -  Parasitic  oscillation  limits  shown  on  one  of  the 
pumping  effectiveness  curves  of  Fig.  19.  The  limitation 
gets  more  severe  as  the  edge  cladding  index  mismatch 
increases.  Also  shown  are  experimental  limits  for  fine- 
ground,  acid-etched  disc  edges. 


amplification  sets  no  definite  upper  limit  on  stored  energy,  but  instead 
makes  pumping  more  and  more  difficult  as  the  energy  density  rises.  Parasitic 
oscillation,  on  the  other  hand,  sets  an  abrupt  upper  limit  to  stored  energy 
at  the  oscillation  threshold. 

Fluorescence  amplification  is  minimized  by  index-matching  the  material 
in  contact  with  the  faces  of  a  laser  disc  to  the  refractive  index  of  the 
disc.  Increasing  the  disc  thickness  causes  only  a  slow  increase  in  the  amount 
of  fluorescence  amplification,  so  thicker  discs  are  of  value  if  pump  uniformity 
and  thermal  distortion  limits  are  not  exceeded. 

Parasitic  oscillation  will  occur  at  very  low  gain  levels  if  the  edge 
cladding  material  has  a  lower  refractive  index  than  the  disc  material.  If 
this  error  is  avoided,  the  level  at  which  oscillation  occurs  will  be 
maximized  if  the  cladding  index  is  as  close  as  possible  to  the  disc  index. 

In  addition,  index-matching  the  faces  raises  the  oscillation  threshold. 
Roughening  of  the  edges  does  not  cause  a  great  increase  of  the  oscillation 
threshold,  and  may  increase  the  average  reflectivity  of  the  edge,  thus 
doing  more  harm  than  good. 

What  limitations  do  fluorescence  amplification  and  parasitic  oscillation 
place  on  glass  laser  amplifier  disc  size?  Figure  20  shows  that  parasitic 
oscillation  is  the  worst  problem.  Unless  stringent  precautions  are  taken, 

oscillation  will  set  in  at  an  across-disc  gain  of  exp(3)  or  less,  assuming  the 

-3 

disc  is  in  air  to  avoid  self- focusing  problems.  With  .5  J  cm  stored,  the 
line-peak  gain  will  be  on  the  order  of  .1  cm  \  This  implies  a  30  cm  disc 
at  most.  This  value  will  vary  with  stored  energy  and  the  amount  of  gain  per 
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unit  stored  energy,  but  values  much  smaller  than  those  quoted  increase 
the  amount  of  energy  lost  to  the  fixed  loss  coefficient  in  the  glass. 

Thus  operation  with  discs  of  dimension  larger  than  about  30  cm  will  be 
very  difficult,  and  operation  in  the  30  cm  region  will  require  great  care 
in  edge- treatment  to  suppress  parasitic  oscillations. 
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